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Abstract 



^^ ! Motivated by recent interest in soft J/^ production in B decays, we investigate B^ 



J/'^D^*' and i]cD^ decays in perturbative QCD. We find that, within that framework, 
these decays are calculable since the heavy cc pair in the final states is created by a hard 
gluon. The branching ratios are estimated to be around 10~^ ~ 10~^, too small to be 
consistent with the data, suggesting that other mechanism(s) contribute to the observed 



o 

Q^' excess of soft J /^ in B^ -^ J/^ + X decays. The possibility of the production of a hybrid 

^^. sdg meson with mass about 2 GeV is briefly entertained. 
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With the advent of the BaBar and Belle B factories, many B decay modes could be studied 
in detail. The rich phenomena of B decays will provide testing grounds for theories of weak 
interactions and hadrons. It is interesting to note that measurements of the inclusive B — *• 
J I'^X spectrum by CLEO [^ and recently by Belle 0, indicate a hump for low J/\E' momentum, 
which kinematically corresponds to J/\& recoiling against a partner as heavy as ~ 2 GeV. 
Brodsky and Navarra [Q suggest that the J /'^ hump may be due to the decay B^ -^ J/'^Ap 
with possible formation oi A — p bound state (an exotic strange baryonium). 

From another view point, Chang and Hou [H proposed as an explanation the existence of 
intrinsic charm in the B meson which decays as B^{dbcc) -^ J/^fD^*^ (and similarly for 77^ 
instead of J/'^)- Thus the intrinsic charm pair transforms into a cc final state while the h 
decays. It is argued that a rate of ~ 10~^ may be possible in this way if the intrinsic charm 
content of B is not much less than 1%. 

We raise here another possibility: The B may decay into a charmonium plus a hybrid, 
B^ -^ J/'^H, where if is a hybrid sdg ^ with Mh ~ 2 GeV 0. Two diagrams that contribute 
to such a process are depicted in Fig. 1. Note that the gluons exchanged in Fig. 1 are soft while 
those in Fig. 2 {i.e. for the conventional B^ -^ J/'^D^*\ see below) are hard, thus enhancing 
the hybrid option as compared to the conventional approach for B^ -^ J/'^D^*\ In addition, 
as shown below, each Feynman diagram in Fig. 2 involves one fermion and one hard gluon 
propagator with average virtuality as large as 10 GeV^. So, we can expect the 5° -^ J /^!H 
decay rate to be 10^ ~ 10^ times larger than B^ -^ J /^!D^*\ although a reliable quantitative 
estimate of the decay rate is very difficult. 

J/* J/* 



B 




(a) 



(b) 



Figure 1: Examples of Feynman diagrams for the production of a hybrid H = sdg in B^ 
J/^H. 



To make such "exotic" suggestions more reliable, one should be convinced that the conven- 



tional picture of heavy mesons indeed leads to tiny numbers in disagreement with experiment. 
To our knowledge, such study is still not available in the literature. In this article we investigate 
these decays within the conventional picture of heavy mesons having the minimal number of 
quarks and using perturbative QCD (pQCD). The applicability of pQCD is justified by the 
large virtuality of the hard gluon which creates a cc pair. As known, in many applications of 
pQCD to B decays |P], say B -^ tttt, the virtuality of the gluon in the hard kernel scales like 
k'l ^ — MbAqcd^; — — 2x GeV^, where x is the momentum fraction carried by the light spec- 
tator quark in the final light meson. However in the processes discussed in this paper (see Fig. 
2), the gluon virtuality scales as k'l > (2?Tic)^. Furthermore, under the common assumption of 
factorization, there are no infrared divergences which cannot be absorbed in wave function, or 
large end point contributions. 
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Figure 2: Feynman diagrams for B^ -^ J/^D*^*' in pQCD. 



We begin our calculation of the decays B^ -^ J/\1'Z)^*-* within the pQCD approach for 
exclusive QCD processes [^ as depicted in Fig. 2, by writing the weak effective Hamiltonian 
Hes for the b -^ cud transitions as 



^eff - -yK^cbV*a 



Ci(/i)rf7^(l -75)mc7'X1 -75)& + C2(/i)c7/.(l -75)wc?7'^(l -75)& , (1) 



where the Wilson coefficients Ci^2(/^) are evaluated to be Ci(mb) = 1.124 and C2{mi,) = —0.273 
at the scale fi = mf, = 4.8 GeV p^ . 

To calculate the amplitudes of the Feynman diagrams in Fig. 2, we take the wave functions 
for 5°, J/^ and D* as follows |n|, g 






4iV, 



4iV. 



it{Mv+ rv)M^)fv, {V = J/^,D*). 



(2) 
(3) 



Now we can write the amplitudes of Fig. 2. as 



I Cf ^ r ,,.,..,. . _ 1 1 



Ma = -ifBf<i>fD*T^as — j^Cijdxdydz(t)B{,x)(t)^,{y)(t)D*{,z) 

xTx[{pB + Mb)757/.(1 - 75) ((1 - x) A- jk + m,)7,] 



1 Cf ^ /-,,,,_,_ , _ 1 1 



D,k^ 



Mb = -?/B/*/D*vra^— -^Ci / dxdydz(t)B{x)(t)^,{y)(t)D*{z 

16 A/^^ J ub 

xTy[{^b + Mb)757.(/; - X /Pb)7m(1 - 75)] 

xTr[/^(Af^+ A)7'^(l - 75) /d-(Mz5.+ /Pd071, (5) 

1 /(^ C \ /■ 11 

Mc = ~ifBhfD*Tras— ( ^Ci + ^C2 J y dxdydz(t)B{x)(t)^,{y)(l)D*{,z) — — 

xTr[(/^B + Mb)757/.(1 - 75)]Tr[/*(M^+ A)7^( A+ /^ + '^c) 

7^^(1-75) itD'{MD^+ ^D*)l% (6) 

-Md = ifBhlD'T^asY^ ( ^C*! + ^^72 J J dxdydz(j)B{x)(f)q,{y)(f)D*{^)^^ 
xTrii^B + Mb)757m(1 - 75)]Tr[/*(M^+ A)7''(l - 75) 
((1 - z) ^n^+ ^) 7. /d*(M,,*+ /P^07l, (7) 

where Cp = | and A^c = 3 is the number of colors. Di{i = a, b, c, d) and /c^ denote the virtuality 
of quark and gluon propagators in Fig. 2.i, which are given by 

Da = -ml + Ml{l-x-y){l-x-z) + {y-z){MUx + y-l)~Ml,{x + z-l)) 

-ze, (8) 

Db = MUx-y){x-z) + {y-z){MUy-x) + Ml,{x~z))-ie, (9) 

D, = -ml + MUl-z) + iMl-Ml,{l-z))z-ie, (10) 

Dd = ^iMUl + 2y-z) + Ml,{l-2y + z) + MU2y^-2yz + z-l))-te, (11) 

e = Mly{y-z)+z{Mly + Ml,{z-y))-ie. (12) 

It may be instructive to evaluate typical virtualities of the propagators involved in pQCD 
calculations. Taking x = 1 — nib/MB, y = 1/2 and z = rric/MB,*, we find 

Da = -20.4 GeV^ Db = 7.2 GeV^ D, = 20.2 GeV^ D^ = 16.4 GeV^ P = 9.9 GeVl (13) 

These values are large enough to justify our pQCD calculation. 



The amplitude for B^ -^ J /'^D* is decomposed as 

A{B -^ J/^{P^)D*{Pd*)) = 4£^.(5 g^, + V Pd*,P^u + iV e^.^^P^P^J, (14) 

where the coefficients S, V and V correspond to s, p and d wave amphtudes, respectively, and 
can be evaluated from Eqs. 4 to 7. The helicity amplitudes are constructed to be 



H, 
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00 



5(M| -Ml- Ml,) + 2PM||p| 



2M,i,Mb. 
H±± = -{S±VMb\p\). (15) 



The branching ratio is given by 



IpI g 



Since the b and c quarks are heavy and their mass is much larger than the typical QCD 
scale Aqcd for a bound state, we can expect that the distribution functions of heavy mesons 
will peak around the points where the heavy quarks are near their mass shell with variance 
^QCT)/fnQ- As an ansatz, the distribution functions are taken as 

0b(x) = 5{x - xb), (p^iy) = 5{y - y^), (pD*iz) = 5{z - zd*), (17) 

with XB = l-j^,y^ = l and zd* = ^. 

To get numerical results, we use V^b = 0.04, fe = 180 MeV, U = 400 MeV, fo* = 
230 MeV, nib = 4.8 GeV, M^o = 5.27 GeV, m, = 1.4 GeV, M^ = 3.1 GeV, M^* = 
2 GeV, a,(2mc) = 0.266. We get 

Br(P° -^ J/^D*) = 6.46 x 10~^ (18) 

and the longitudinal polarization fraction is 

Pl = y = 0-398- (19) 

Since the amplitudes are highly suppressed by the large virtualities of the propagators as shown 
in Eqs. 8 to 13, the smallness of Br(P° -^ J/'^D*) is understandable. To illustrate the stability 
of our results, we plot in Fig. 3 Br(P*' -^ J/'^D*) versus xb-, «-e., the peak point of (J)b{x)- 




Figure 3: Br(i?° -^ J j'^D*) vs. xb, the peak point of 0b(x). 



From Fig. 3, we can see that the rate is rather stable against changes of the parameter xb- 
Due to relativistic effects, the distribution functions should have variances of O^Aqcb/i^q)- 
To show the effects of the variances, we take 



^b[x 



^D'[X 



H[x 



Nbx{1 — x) exp 
Nd^x^I — x) exp 
A^*x(l — x) exp 



f Mb \\ ., 

- -T7 [x - Xb) 

\Mn - nihJ 



Mb -nib. 



M, 



D* 



nir 



Mvp 



ix - Xd* 



X — 



(20) 
(21) 
(22) 



where Nb, No* and A'^^ are normalization constants to make J dx(f){x) = 1. To model the 
distribution functions, we take the mass difference between the heavy meson and its heavy 
constituent (s) as shape parameter. These distribution functions follow the consensus that the 
smaller the mass difference the sharper the distribution functions. Using these distribution 
functions, we obtain 



Br(5° -^ J/^D*) = 8.50 x 10" 



Pr 



T 



0.395. 



(23) 



Since B^ decays with three charm quarks in its final states, it could be taken as a probe 
of strong interactions, especially hadron dynamics. We extend our calculations to B^ — >■ 
J/'^D, TjcD and t]cD* decays. The amplitudes for these decays can be obtained through the 
following replacements in Eqs. 3 to 7 

i 









Vc 



AN, 



AN 



■75(/'d +MD)Mx)fD, 
■75(Ac +Mr,J<l)r,MU- 



(24) 



Using fi) = 200 MeV, /^^ = 335 MeV, the branching ratios are estimated to be 

Br(fi° -^ J/^D) = 7.28 x 10-^ 
Br(fi° -^ r]cD*) = 1.39 x 10"^ 
Br(fi° -^ 7]^D) = 1.52 X 10"^. (25) 

In summary, we have studied the decays B^ -^ J/'^{ric)D^*^ within the conventional theoret- 
ical framework. The branching ratios of these decays are estimated to be around 10~^ ~ 10^^. 
B^ decays to J/^!D^*' can not account for the excess for slow J/\l/ as indicated by the CLEO 
measurement of the J/\l/ momentum spectrum in B inclusive decays. Experimentally, inclusive 
decays of B mesons to charmonium could be well studied at BaBar and Belle, and it is impor- 
tant to confirm whether the slow J/\E^ hump exists with refined measurements. If the excess 
persists, it would be hard to explain the phenomena within the conventional theoretical frame- 
work for hadron dynamics. As shown in here, our numerical results are rather stable under the 
change of parameters. // these exclusive decays were observed to be abnormally large, say, of 
order 10~^ ~ 10~^, it would challenge the conventional theoretical framework and bring forth 
new interesting QCD phenomena, like the scenarios discussed in Ref. 0, ^ or the possibility 
raised here, of the formation of a ~ 2 GeV sdg hybrid state H through B^ -^ J/'^H. Finally 
let us note, that multibody final states such as J/^D^*^ + nvr, where n^ax = 1? 2 for D*, D, 
respectively, being on the edge of phase space, are expected to be even smaller than those with 
n = 0. 
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